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Abstract 



Lower bounds on the magnitude of the spectrum of the Hermitian Wilson- 
Dirac operator H{m) have previously been derived for < m < 2 when the 
lattice gauge field satisfies a certain smoothness condition. In this paper lower 
bounds are derived for 2p — 2<m<2p for general p = 1,2, ... ,d where d is the 
spacetime dimension. The bounds can alternatively be viewed as localisation 
bounds on the real spectrum of the usual Wilson-Dirac operator. They are 
needed for the rigorous evaluation of the classical continuum limit of the axial 
anomaly and index of the overlap Dirac operator at general values of m, and 
provide information on the topological phase structure of overlap fermions. 
They are also useful for understanding the instanton size-dependence of the 
real spectrum of the Wilson-Dirac operator in an instanton background. 

1 Introduction 

It is well-known from numerical studies (see, e.g., [1, 2]) that in smooth gauge back- 
grounds in d dimensions the real eigenvalues of the Wilson-Dirac operator are localised 
around the values 0,2,4, ... ,2d (in units of inverse lattice spacing, and with Wilson 
parameter r = 1). In this paper we give an analytic derivation of this numerical 
observation. Our smoothness condition is the "admissibility condition" of [3, 4]: 



1-U{p)\\ < e 



Vplaquette p. 
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Since the plaquette variable has the expansion U{p) — 1 — a^F^,^(x) + 0(a^) in powers 
of lattice spacing a, (1.1) can be regarded as an approximate smoothness requirement 
on the curvature of the lattice gauge field. If U is the lattice transcript of a smooth 
continuum gauge field then (1.1) is automatically satisfied for any e > when the 
lattice is sufficiently fine. 

In fermionic definitions of the topological charge of lattice gauge fields the low- 
lying real eigenmodes of the Wilson-Dirac operator Du, are interpreted as would-be 
zero-modes, while the other real eigenmodes are interpreted as would-be doubler 
modes. This interpretation rehes on the real eigenvalues being locahsed as described 
above, which is not the case in general for arbitrary rough gauge fields. The local- 
isation result on the real spectrum of derived in this paper provides a specific 
analytic criterion under which the localisation is guaranteed. It is also of interest in 
connection with the overlap fermion formulation on the lattice [5, 6]. This is because 
a real eigenmode for the Wilson-Dirac operator is equivalent to a zero-mode for the 
Hermitian Wilson-Dirac operator with negative mass parameter: 
in 

D^ip = — -0 <^ H(mU = 75(aD^ - m)ip = (1.2) 
a 

Localisation of real eigenvalues of D^, around 0, 2, . . . , 2ci (in units of 1/a) is therefore 
equivalent to absence of zero- modes for H{m), i.e. to existence of non-zero lower 
bounds on \H{m)\, when m is away from these values. This implies a topological 
phase structure for the overlap Dirac operator [6] Dov = ^(l + 75-ff(m)/|if(m)|), since 
the index of D^v (a well-defined integer) is locally independent of m but can jump at 
the values for which H{m) has zero-mode(s). The topological phase structure of 
has previously been studied in ref.'s [7, 8]. The bounds derived in this paper lead 
to analytic information on the topological phases which complements the numerical 
results of those papers. 

Furthermore, a non-zero lower bound on \H{m)\ allows the locality of the overlap 
Dirac operator, and its smooth dependence gauge field, to be analytically established 
[3] (see also [9]). The general bounds derived in this paper allow the unnatural 
restriction < m < 2 on the results of [3] to be removed. These bounds are also 
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required for the rigorous evaluation of the classical continuum limits of the axial 
anomaly and index of the overlap Dirac operator [10, 11].^ As a final application we 
will discuss qualitative implications of the bounds for the instanton size-dependence 
of the real spectrum of the Wilson-Dirac operator in an instanton background. 

The paper is organised as follows. In Section 2 the previously derived lower 
bounds on \H{m)\ are summarised and the new general bounds are formulated. The 
new bounds are derived in Section 3. The derivation is rather technical and not so 
illuminating, so in Section 4 we supplement it with with a heuristic argument which 
provides a clearer intuitive understanding of why the bounds exist. The heuristic 
considerations are further developed to give an analytic explanation of properties of 
the spectral flow of H{m) previously observed in numerical studies. In Section 5 
the above-mentioned applications of the bounds are discussed, and the results of the 
paper are summarised in Section 6. A generalisation of the bounds from the standard 
case of Wilson parameter r = 1 to general values r > is given in an appendix. 

2 Summary of previous bounds and formulation of the new 
bounds 

For m < and m>2d ((i=spacetime dimension) it is well-known that \H{m)\ > \m\ 
and \H{m)\ > m—2d, respectively; see, e.g., [5, 17, 2]. (We review these bounds and 
generalise them to arbitrary values of the Wilson parameter r in the Appendix.) By 
(1.2) this implies that the real eigenvalues of D^, (in units of 1/a) must lie in the 
interval [0,2d]. 

In [3, 4] lower bounds of the form 

\H{m)\ > Vl-cie- |l-m| (2.1) 

were derived when the lattice gauge field satisfies the smoothness condition (1.1). The 

currently sharpest bound has ci = 6(2+\/2) 20.5 in 4 dimensions [4] and generalises 
^ Other evaluations of the classical continuum limit of the axial anomaly (less rigorous, and not 
using a lower bound on \H{m)\) have been given in [12, 13, 14, 15]. 
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to Ci — {2 + ^/2)d{d — l)/2 in d dimensions. Clearly (2.1) can only be a nontrivial 
lower bound if e < 1/ci and |1 — m| < a/1 — Cie. The latter implies < m < 2. Lower 
bounds on \H{m) \ in the "doubler regions" 2 < m < 4, 4 < m < 6, ... 2d— 2 < m < 2d 
have so far been missing. 

Note that by (1.2) existence of a nontrivial lower bound on \H(m)\ for |1 — m| < 
a/1 — Ci e is equivalent to the Wilson-Dirac operator having no real eigenvalues 
in the open interval ]1 — a/1 — Ci e , 1 + ^/l — Ci e [. To extend this to a general 
localisation result for the real eigenvalues of D^j existence of lower bounds on \H{m)\ 
ior k — 1 < m < k+1 , k — 1,3, ... , 2d—l needs to be established. 

Our aim in this paper is to generalise (2.1) to bounds of the following form: 

\H{m)\ > Vl-Cfce- \k-m\ , k = 1,3,5, ... ,2d - 1 (2.2) 

For given m & ]k — l , k + l[ this lower bound is nontrivial when e in the smoothness 
condition (1.1) is chosen such that e < (1 — (A; — m)^)/cfc. On the other hand, if we 
only require e < l/c^ for all k then the bound is nontrivial for all values of m except 
those lying in one of the following intervals: 

[0,1- vT^^ ] , 

[k + VI- Cke , k + 2- ^1- Ck+2 e ] , k^l,3, . . . ,2d-l , 
[2d-l + ^1 - C2d-i e , 2d ] (2.3) 

illustrated in Fig. 1 below. In this case the real eigenvalues of (in units of 1/a) 
must lie in these intervals. Clearly when e is small these intervals are localised around 
the values 0, 2, 4, . . . , 2ci. More specifically, we see that the real eigenvalues of Dyj are 
guaranteed to lie in the intervals [0,S], [2p — S,2p + S] {p—1,2, . . . ,d — l), [2d — S,2d] 
when e < (1 — (1 — 5)^)/cfc for all A; = l, 3, . . . , 2d — 1. This is the advertised localisation 
result on the real spectrum of Dyj. Explicit values for the c^'s will be determined in 
the next section. 
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1+8, 8 



Figure 1: Illustration of the intervals (2.3) (with Sk = ^yl — e) in the dimension 
d — A case. The bounds (2.2) imply that the real eigenvalues of Dyj he in these 
intervals. 

3 Derivation of the bounds 

The Wilson-Dirac operator D^, with general Wilson parameter r is given by 

aD^{T)^Y.l'\iT+^^-T-,)+r{l-\{T^, + T_,)) (3.1) 

where T^^ are the forward/backward parallel transporters ((T+^)a;j, = Un{x)Sx,y-fi , 
(2±(u)* = {T±ij)~^ — T^fj). Dyj is an operator on the lattice spinor fields living on a 
hypercubic lattice on an even d-dimensional Euclidean spacetime and taking values 
in some unitary representation of the (unspecified) gauge group. The spacetime may 
be either the infinite volume R'^ or finite volume rf-torus T'^. In the former case the 
(completion of the) space of spinor fields is an infinite-dimensional Hilbert space, 
while in the latter case it is simply a finite-dimensional complex vector space with 
inner product. In the following || • || denotes operator norm. Clearly ||T±^|| = 1, so 
Dyj is bounded. A well-known, important consequence of (1.1) is 

II [T±,,T^u] II < e , II [T±^,T^.] II < e (3.2) 
It is useful to define the hermitian operators 

These have bounds — 1<5'^<1,— l<C^<l,0<i?^<2 and satisfy (in any 
gauge background) the following identities: 

[S^,C,]^0 , [S,,R,]^0 (3.4) 

Sl + Cl = l , Sl = R,{2-R,) (3.5) 
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The Wilson-Dirac operator can then be written as 



Dn.{r)^-Yl i^S^ + rR^ (3.6) 



For later use we also note the relations 



R, = IK'^, 2-i?, = l(2 + V^r(2 + V,) (3.7) 



where V^^ = - 1. 

The Hermitian Wilson-Dirac operator (normalised by 1/a) is given by 

Hi^rn, r) = ^^{aDwir) — rm) (3.8) 



\H{m,r)\ = y H{m,ry is defined via spectral theory. In the following we set r = 1 
and consider H{m) = H{m, 1); the case of general r is dealt with in the Appendix. 

To derive the desired bounds (2.2) it suffices to show the following: 

H{kf > 1-Cfee , A; = l,3,...,2d-1 (3.9) 

Indeed, the eigenvalues A(m) of H{m) satisfy | < 1 [16, 4], implying \H{rn')\ > 
\H{m) \ — |m — m'\ (an alternative derivation of this had also been given in the first 
paper of [3]), and this together with (3.9) implies the bounds (2.2). 

To derive bounds of the form (3.9) we use (3.5)-(3.7) to express H{mY as follows: 

H{mY = {aD^ — m){aD^ — m) 

= 1 + X{m) + E' (3.10) 



where 



and 



X{m) = j:Sl + {-m + j:R,r-l (3.11) 
= ^i?^i?,-2(m-l)^i?/, + m'-l (3.12) 



E' = ^rYl[s,,s.]+ir[s,, a] (3.13) 
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Using (3.2) and triangle inequalities a bound on E' of the form 

1 1^1 1 < c'e (3.14) 

can be obtained. The value for d obtained in [4] in the 4-dimensional case is d — 
6(1 + \/2) ~ 14.5 and generalises to c' = (1 + \/2)d{d — l)/2 in d dimensions. 

To complete the derivation of (3.9) we need to show that x(A;) can be written in 
the form 

X{k) ^ P{k) + E{k) , P{k) > , \\E{k)\\ < 4e for A; = l,3,...,2ci-1 

(3.15) 

It then follows from (3.10) that (3.9) is satisfied with Ck — c'l + c'. 

It is easy to derive a decomposition and bound (3.15) in the k — 1 case [3, 4]. In 
this case (3.12) reduces to 

X(1) = E^M^- (3-16) 

Using (3.7) one finds R^R, = | V^V^i?, = | V^i^^V^ + ^ V;[V^ , R,] ^ P^, + E^, 
where P^^ = IW^R^V ^ > and \\E^^\\ < e, leading to x(l) = -P(l) + E{1) with 
P(l) > and ||£'(1)|| < d{d — l)e in d dimensions [3]. A more subtle decomposition 
x(l) = -P(l) + was derived in [4] for which ||£'(1)|| < ^d{d — l)e. In this way 
the k ^ 1 bound (2.1) was obtained with ci = c'/ + c' = 6 + 6(1 + s/O.) fa 20.5 in 4 
dimensions [4]. 

Our goal now is to derive a decomposition and bound (3.15) for x{k) in the case 
of general k = 1,3, . . . ,2d — 1. Setting 

4°) = 2-i?^ and (3.17) 

we begin by noting the identity 

X{m) = x{m) + x{fn)rev (3.18) 

where 

X(m) = -1^ E ((m - 2(gi + . . . + q,)f - l)Rt^ R^^ ■ ■ ■ R^^^ (3.19) 

9i,...,ijd=0,l 
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X{m)rev is defined by replacing R'f'^Ei^^^ ■ ■ ■ B^i"'^ by B^i"'^ ■ ■ ■ R^^'^R^^'^ in (3.19). The 
key feature of this expression is that, unhke the original expression (3.12), it is a sum 
of monomials in the positive operators i?^ and 2 — Riy (recall < i?^ < 2) with positive 
coefficients when m is an odd integer (in particular when m = /c= 1, 3, . . . , 2d—l). As 
we will see shortly, this provides for a decomposition x{k) — P{k) + E{k) of the form 
required in (3.15). 

To derive (3.18), consider the expansion of xiv^) powers of the i^^^'s: 
^ = aQ + ai'^R^ + . . . + ap ^ R^^- ■ ■ R^,^ + . . . + adRiR2- ■ ■ Rd (3.20) 

H iJ,i<—<tJ,p 

The expansion of x{''^)rev is identical except that the ordering of the R^^s is reversed. 
In light of (3.12), to derive (3.18) it suffices to show that 

ao = -^{j^^ ~ 1) ) q;i = — (m — 1) , q;2 = 1 and ctp = for p > 3 (3.21) 

Let us focus on the term of order p in (3.20). It gets contributions from the terms in 
(3.19) with ?i + . . . + ?d < p. The terms with gi + . . . + g-j^ = s are 

2^)2-1) (2-i?i)---(2-i?.,_i)i?.,(2-i?,,+i)-.. 

(2 - i?.,_i)i?.,(2 - it:.^+i) • • • (2 - i?,) (3.22) 



For s <p the contribution of this to the ctp-term in (3.20) is 
^((m-2.)2-l)(-l)^-2'^-^[^] 



(3.23) 



/ii<---<//p 

(the binomial coefficient ^ appears because it is the number of ways to pick s distinct 
elements from a set of p elements) . It follows that 

1 P 



Y{{m-2sf-l) 



s=0 



.ly- 



(3.24) 



From this we find cto — 2(^^ — ^) ■, oti — — (m — 1) and a2 — 1 as claimed in (3.21). 
In the p > 3 case we calculate 



^ ^(4^2 - 4ms + - 1) M (-1)^' 



s=0 



+(^'-i)EH (-ir 



p-i 



s=0 s=0 

P 



s=0 



(3.25) 



(each sum vanishes since X^s^o (~1)^ ^ '* = (! — 1)^ ^ etc.). This completes the 

derivation of (3.21), thereby estabhshing (3.18). 

We now show how (3.18)-(3.19) leads to a decomposition x{^) = Pi}') + -^(^") 
of the form (3.15). The operator product • • ■ R^^'^'' in (3.19) decomposes into a 
positive piece and a piece involving commutators as follows. Set 

V(f) -2 + V^ = r+^ + l and V^^^ - = - 1 (3.26) 

then ||Vj«)|| < 2 for g = 0, 1 and, by (3.7) and (3.17), 

R^^ = \iyfr^f- (3.27) 

Using this and commutator relations [O, Oi • • • Op] = YZ=i Oi • • • Os_i[0, Os]Os+i ■ ■ ■ Op 
we obtain 

^(gi) . . . ^fe) _ p(9i,...,9,) ^ (3 28) 

with 

p(gi,...,9d) ^ J_(y{?ci)...y(9i))*yfe)...y{?i) (g^Q) 

2*^ 

X ( E 4''^ • • • ^?-r'^ [v^^ ' ] • • • ^^Vi-r'^ • • • vS''^^ 

^ s=p ' 

(3.30) 



Clearly p(«i'-'«d) > 0. Furthermore, the bounds ||Vjf'')|| , \\R^^^W < 2 and, by (3.2) 
II [^If^^ -Ri?"'* ] II ^ ^) together with triangle inequalities, lead to the bound 

p(9i,...,9d)|| < (3.31) 
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where 

d-l d 

c=J2Yl 2'^"^"^ 2f-^ = 2"^-^ {d -l){d + 2) (3.32) 

p=l s=p 

The reversed product R^^'^^ ■ ■ ■ b!i^^ has an analogous decomposition P^^^'-''^'^^+E^ll'-''^''^ 
with identical bounds. Consequently, by (3.18)-(3.19) we get a decomposition 

X(m) = P(m) + E{m) (3.33) 

where P(m) and E{m) are given by (3.19) with R^^^ ■ ■ ■ R^^"^^ replaced by p(9i'---'9d) -)_ 
P^iv'"''^'^^ and £;(9i. -.9d) _|_ £!(9i. -.9d)^ respectively. The coefficient in the summand in 
(3.19) is > when m is an odd integer, hence P{k) > for odd k and in particular 
for = 1, 3, . . . , 2d — 1 as required in (3.15). Furthermore, from (3.31)-(3.32) we get 
a bound 

\\Em < 4e (3.34) 

with 

^ ^ qi,...,qd.=0,l ^ 
^ p=0 

= 2'^-^{d-\){d^2){{k-df -\^d) {d>2) (3.35) 

Thus we have established the existence of a decomposition and bound (3.15) for xik) 
for general k = l,3,...,2d— 1. By our previous discussion this implies the existence 
of the desired bounds (2.2). We remark that (3.35) is invariant under k ^ 2d — k. 
This is as expected in light of the well-known fact that a lower bound on \H{m)\ is 
also a lower bound on \H(2d — m)\ (see the Appendix). 

The bound (3.34)-(3.35) is rather weak. For example, in the c? = 4 case it is 

cl{d = 4) = 36{{k - 4)^ - 1) + 144 , (3.36) 

giving in the k — 1 case c'( — 432, which is much larger than the values c'/ = 12 and 

c" = 6 obtained in [3] and [4], respectively. Note however that for the applications 
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discussed in this paper it suffices simply to sfiow tfie existence of bounds of tfie form 
(2.2) without necessarily finding sharp ones. The largeness of dj, in the above bound 
is due to the large number of terms in the expression (3.19) for x(m). In practice it 
is often possible to simplify this expression such that a sharper bound (i.e. smaller 
dl) can be derived. We discuss this in the d — A case in the following. 

In the remainder of this section we specialise to dimension d = 4 and consider x(/c) 
for k — 1,3,5,7. We wish to simplify the expression (3.18)-(3.19) for x(^) order 
to get bounds with smaller c'^. In order to have a decomposition x(/c) — P{k) + E{k) 
the simplified expression must continue to be a sum of monomials in the positive 
operators R^, (2 — Ry) with positive coefficients. In the k — 1 case (3.18) simplifies 
to x(l) = Hn^u R^J.Ru (recall (3.16)) from which the previously discussed bounds 
with c" = 12 [3] and c'/ = 6 [4] can be derived. In the k = 7 case x(7) reduces to 
E^<.(2 - R^){2 - R,), leading to 

x(7) = E(2-^M)(2-i?.) (3.37) 

Arguments analogous to the ones in [3], and [4], lead to bounds with Cy = c'/ = 12, 
and Cy = c'/ = 6, respectively. Turning now to the A; = 3 case, (3.19) gives 

X(3) = i(2-7?i)(2-i?2)(2-i?3)(2-i?4) + |i?ii?2i?3i?4 

+ I ((2 — R\)R2R^Ra + -Ri(2 — R2)RiRi + RiR2{'2 — R^)Ri 

+RiR2R3{2 - i?4)) (3.38) 

In this case there does not appear to be a major simplification with the required 
properties. In fact it is quite easy to show that x(3) cannot be written as a sum of 
monomials of order < 3 in i?^ , (2 — Ri,) with positive coefficients (we leave this as 
an exercise for the reader). Minor simplifications are possible though, for example: 

X(3) = i((2-i?i)(2-i?2)(2-i?3)(2-i?4) + i?li?2i?3(2-i?4)) 

+ |(i?i-R2-R4 + R1R3R4 + -R2-R3-R4) (3.39) 

x{^)rev simphfies analogously. Estimates of the kind used to derive (3.31)-(3.32) 

show that the decomposition P(3) +£'(3) of the resulting expression for x(3) satisfies 
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ll-E'lS)!! < Cg'e with Cg = 42. This is considerably smaller than the value Cg = 144 
provided by (3.36). It is plausible that a bound with even smaller Cg can be derived, 
e.g. by an extension of the arguments of [4], but we will not pursue this here. Finally, 
in the k = 5 case analogous arguments lead (as expected) to a bound with Cg = Cg = 42 
(we omit the details). 

4 Heuristic considerations 

In this section we present a heuristic argument which provides a clearer intuitive 
understanding of why bounds of the form derived in the previous section should 
hold. We go on to heuristically derive certain properties of the spectral flow of H{m) 
previously observed in numerical studies (e.g. [1, 17]). 

Consider a "near zero- mode" for H{m): 

H{mf^p fti (4.1) 

If e in the smoothness condition (1.1) is small then E' (3.10), and (4.1) becomes 

(recall C^ = l- R^) 

iT.Sl + {-m + E(l - C,)f (4.2) 

Since > it follows that ~ for // = 1, . . . , 2d and consequently, by (3.5), 

Cltl^^il- S^)^ ^ -0 which implies that C^ip ^ (-1)^'^^/^ for j^^O or 1. Then (4.2) 
reduces to 

^ ( -m + $^(l-C^))'V' 

~ ( -^ + IZ(1 - (-iyO)^V' (4.3) 

which implies 

m ^ ^(1 - i-iy^) . (4.4) 

Thus we see heuristically that when e is small the only values of m for which H{m) 
can have "near zero-modes" are m = 0, 2, 4, . . . , 2(i. This makes plausible the result 
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of the previous section, namely that when m is away from these values a nonzero 
lower bound on \H{m)\ should exist. 

In fact the above heuristic approach can be further developed to get an alternative 
rigorous derivation of the bounds (2.2) [18]. However, the argument is technically 
more complicated than the one in §3, and does not lead to sharper bounds, so we do 

not present it here. 

We now proceed to study the spectral flow of H{m). For this it is useful to 
introduce the operators defined by^ 

(Vxy = 757^(-l)"''<^x, (n^ = xja e Z) (4.5) 

These have the following properties: T^ — —l, 7^7^ = lor ji ^ v , [7^ , ] = 

, C^,% = -%C^, , [7^ , ] = for // 7^ i/, and 7^75 = -75 T;^^. Using these we find 

H{m)% = -%{H{m) + 275(7^) (no sum over //) (4.6) 

By (3.2) [H{m) , ] when e is small, so the eigenspaces of H{m) can be decom- 
posed into approximate eigenspaces for the C^^'s. I.e. for eigenvectors ■0(m) of H{m) 
with H{m)ip = \{m)ijj{m) we can assume C^'0(m) ~ c^ij){m). The eigenvalues are 
independent of m since is independent of m and has discrete spectrum. Then, 
by(4.6), 

H{m)Tnip{m) = —TnH{m — 2c,j) ip{m) (no sum over /i) (4.7) 

Set ipf_i{rn) := T^'il){m — 2c^). It follows from (4.7) that ipfj,{m) is an approximate 
eigenvector for H{m) with eigenvalue ~ — A(m — 2c^). Similarly we find 

^M^Mi-M.M ^ (-irA(m-2(c^, + ...c^J)^^,..,^^(m) (4.8) 

where ^^i...^p(m) := 7^^ ■ ■ ■ 7^^ ^(m-2(c^^ + . . .+c^p)), when the //^-'s are all mutually 
distinct. 

^These have proved useful in previous lattice fermion contexts; see, e.g., [13] and the ref.'s therein. 
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Now, if A(m) crosses zero near m — then by our previous argument (recall (4.4)) 
^ {-ly^ with E;u(l - (-1)-^^) ~ 0, i.e. ~ 1 for all /i , and (4.8) becomes 

^(^)V'^i-^.M ~ (-l)^A(m-2p)V'^,...^^(m) (4.9) 

i.e. ipn^...ij,^{m) is an approximate eigenvector for H{m) whose approximate eigenvalue 
A(— l)^A(m — 2p) crosses zero near m — 2p. Furthermore, the sign of the crossing 
is (—1)^ relative to the sign of the crossing of zero by A(m) near m = 0. We note 
the following: (i) If {/xi, . . . , /Xp} 7^ {i^i, . ..,i^p} then V'mi-MfI"^) ip^^...^,^{m) are 
approximately orthogonal since they are approximate eigenvectors for the C^'s with 
different eigenvalues, (ii) ■^^^...^^(m) is unchanged up to a sign under a change of 
ordering of the //j's (since 7^7^ = for 1^ ^ Hence we can assume that the 

jUj s are ordered so that fii < . . . < ^p. (iii) If -(/^(m) is an eigenvector for H{m) whose 
eigenvalue A(m) crosses zero at some value mo then by (4.4) mo ~ Z)ju(l ~ i^^Y'') 
where Cju'0(m) ~ (— l)-''"0(m) for — or 1. Any such eigenvector arises in the 
way described above, i.e. ip{m) — i^f^i-npiiTi) — 7^^ • • ■Tn^%l){m — 2p). Indeed, we 
set -0(771) = (— l)^7^j ■ ■ ■ Tp_^ip{m + 2p) with the /x/s being the /x's for which — 1. 
Then by (4.8) ilj{m) is an approximate eigenvector for H{m) whose eigenvalue \{m) 
is ~ at some value of 777 near zero. (To see this, recall = —1.) Thus we have 
heuristically established the following: The eigenvectors of H{m) whose eigenvalues 
cross zero at some value of m can be naturally grouped into sets of 2'^ elements: 
One of the eigenvectors, ip{m), has eigenvalue X{m) crossing zero near m — with 
crossing sign ±. There are d eigenvectors, ipdlm), with eigenvalues crossing zero near 
777 = 2 with sign ^, and more generally d\/{p\{d — p)!) eigenvectors '^^^...^^(7^7) , ni < 
... < /ip , with eigenvalues crossing zero near m — 2p with crossing sign ^(—1)^ for 
p — 1,2, . . . ,d. This is precisely the spectral flow property of H{m) found in numerical 
studies in 2 and 4 dimensions [1, 17]. An illustration of the spectral flow associated 
with one such family in the (i = 4 case is given in Fig. 2. The Hermitian Wilson- 
Dirac operator in any gauge background U has the well-known property H{U,m) — 
—H{—U, 2d — m), so that H{m) — OH{U, 2d — m)0~^ for a certain unitary operator 
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Figure 2: Illustration of the spectral flow associated with a typical eigenvector family 
of H{m) of the kind discussed in the text in dimension (i = 4: A crossing near m = 
with +sign is associated with 4 crossings near m — 2 with —sign, 6 crossings near 
m = 4 with +sign, 4 crossings near m = 6 with —sign and 1 crossing near m = 8 with 
+sign. 

O (see, e.g., [2, 17]).^ Hence if A(m) is an eigenvalue for H{m) then — A(m) is an 
eigenvalue for H(2d — m). This property must be manifested in the eigenvalues of 
the family of eigenvectors of H[m) discussed above, and is also illustrated in Fig. 
2. Combining this spectral property of H{m) with the fact that the index of the 
overlap Dirac operator equals —1/2 times the spectral asymmetry of H{m) [5, 6], an 
immediate consequence is the relation index{Dov{m)) — —index{Dov{'2d — m)) which 
was emphasised in [7]. 

5 Applications of the bounds 

We have already seen in §2 how the general bounds (2.2) lead to a localisation re- 
sult on the real spectrum of the Wilson-Dirac operator, thus providing an analytic 
understanding of the numerical results on the real spectrum in "smooth" gauge back- 
grounds. In this section we discuss applications of the bounds to overlap fermions [6]. 
The general bounds allow analytic results on the overlap Dirac operator Dgv which 
^For O to exist the number of lattice sites along each edge of T'' must be even. 
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were previously derived for the < m < 2 to be extended to the general m case 
(m 7^ 0,2,4, .. . ,2d). Although < m < 2 is the physically relevant case (i.e the 
case where Dov is free from spurious fermion species) this restriction appears quite 
unnatural and it is of some theoretical interest to know the properties of Dgv in the 
regions 2p < m < 2p + 2 , p—1,2, . . . ,d where the extra fermion species are present. 

(i) Locality and smooth gauge field- dependence of the overlap Dirac operator. 

With the bounds (2.2) the arguments of ref.[3] for the locality of Dov and its 
smooth dependence on the lattice gauge field carry over unchanged from the < 
m < 2 case to the k — l<m<k+l case {k = 1,3, . . . ,2d — 1) after choosing 
e < {1 — {k — mY)/ck so that the lower bound on \H{m)\ is greater than zero. The 
size of the exponential decay constant in the locality bound for Dov depends on the 
size of Ck , but for the existence of the locality bound it is enough to know that (2.2) 
holds for a specific value of which is independent of the lattice gauge field. 

(ii) Evaluation of the classical continuum limit of the axial anomaly and index of the 
overlap Dirac operator. 

The rigorous evaluation of the classical continuum limit of the axial anomaly^ 
and index of the overlap Dirac operator at general values of m requires existence 
of a nontrivial lower bound on |i/(m)| when the lattice is sufficiently fine [10, 11]. 
We claimed in [10, 11] that such bounds exist and promised to provide them in a 
forthcoming paper. The present paper delivers on that promise. Again, the actual 
values of the c^'s do not matter: The lattice transcript of a smooth continuum gauge 
field automatically satisfies the smoothness condition (1.1) for any e > when the 
lattice is sufficiently fine (see [11] for the rigorous justification of this point), so all 
that matters for the classical continuum limit calculations is that the bounds hold for 
some choice of c^'s which are independent of the gauge field and lattice spacing. 

*It may sound contradictory to speak of the "classical" continuum limit of a purely quantum 
quantity such as the axial anomaly, so let us explain the meaning: "Classical" refers to the fact that 
one considers the a — > limit of the axial anomaly with the lattice gauge field given by the lattice 
transcript of a smooth continuum gauge field. See [10, 11] for the details. 
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(iii) Topological phase structure of the overlap Dirac operator. 

In the finite volume ci-torus case the index of Dgv = ^(1 + j^H {m) / \ H {m)\) is a 
well-defined integer; it is locally constant in m but may jump at the values at which 
H{m) has zero-mode(s) . Thus Dgv has different topological phases and the value of 
m should be chosen so that D^^ is in the "correct" phase. This issue has previously 
been studied both analytically and numerically in [7] and numerically in [8]. However, 
the analytic arguments in [7] arc problematic since they involve treating topologically 
nontrivial fields as perturbations of the trivial gauge field U—1. On the other hand, 
the bounds (2.2) provide rigorous nonperturbative insight into the topological phase 
structure when the lattice gauge fields are required to satisfy the smoothness condition 
(1.1) with e < 1/cfe for all /c = 1,3, . . . ,2d — 1 : they imply that there are distinct 
topological phases for D^y , with each phase characterised by m being in one of the 
open intervals ]k — — Ck€ , k + — Ck€ [. The result of [11] states that for SU(A'") 
gauge fields on the d-torus {d — 2n, n > 1), or U(l) gauge fields on the 2-torus, 
index (Dov) coincides with the index of the continuum Dirac operator in the classical 
continuum limit provided < m < 2.^ This indicates that the "proper" topological 
phase for Dov is the one where m is in the interval ]1 — \/l — Ci e , 1 + -s/l — Cie [. We 
denote index{Dov) by Q when D^v is in this phase. A complete description of the 
topological phases for Dov when the smoothness condition (1.1) is imposed is now as 
follows. For m < , is in a topologically trivial phase (i.e. index{Dov) = in any 
gauge background) [5]. For < m < l — — Cye , D^v is not in a distinct topological 
phase: index{Doy) can be any value from to Q depending on the background gauge 
field. In a given gauge background, as m is increased from to 1 — y/1 — ci e , the total 
spectral fiow of H{m) is Q. This is due to the well-known fact that at each crossing of 
zero by an eigenvalue of H{m) the index of Dov changes by =f1 depending on the sign of 
the crossing. For I — ^1 — ci e <m < l + — ci e , Dov is in the "proper" topological 
phase where index{Dov) — Q- For 1 -|- ^/l — ci e < m < 3 — ^/T^^C3e , Dov is no 
longer in a distinct topological phase and the spectral fiow of H{m) as m increase 
^This was shown in [11] in the case of the 4-torus, but the argument generalises straightforwardly 
to the general rf=2n-torus. 
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through this region is —dQ. For 3 — — C3 e < m < 3 + -\/l — C3 e , Dot, is in another 
distinct topological phase with index{Dov) = (1 — d)Q. The pattern continues as m 
increases: For k—y/1 — Ckt < m < k+y/1 — c^e , Dqv is in a distinct topological phase 
with index{Dov) — ( X]p='o^''^^(— 1)^ p ) Then, as m increases from k + ^1 — e 
to A; + 2 — y/1 — Ck+2 e , -Dot; is no longer in a distinct topological phase and the 
spectral flow of H(m) through this region is (— ^f,_^^y^ Q. Finally, after m 
has increased to 2d we have index{Dov) = ( EJ=o(-1)*' [p]) Q = - '^Y Q ^ ^ and 
Dgy is back in a topologically trivial phase which it remains in for all m > 2d. 

The above description of the topological phase structure of Dg^ is compatible 
with the results of previous numerical studies in 2 and 4 dimensions [7, 8]. To put 
the above analytic argument on a completely rigorous footing a rigorous derivation 
of the heuristic result of §4 on the spectral flow of H{m) is required. This remains 
as a problem for future work. We note however that further evidence for the validity 
of this description comes from the result of [11], which states that in the classical 
continuum limit 

(fc-l)/2 

index{Dov) = ( ^ (-1)^ [^j ) Q for A;-l< m < A;+l (A; = l, 3, . . . , 2d-l) 

= for m < and m >2d 

where Q is the index of the continuum Dirac operator. 

A generalisation of the overlap Dirac operator has been presented in [19] and it 
would also be of interest to estabhsh the topological phase structure of this operator. 
For this, bounds on the generalised Hermitian Wilson-Dirac operator for general m are 
needed. A bound has already been derived in [20] for < m < 2 by a generalisation 
of the argument of [4]. It is plausible that bounds for general m can be derived by 
a generalisation of the argument in the present paper. We leave this as a potential 
topic for future work. 

(iv) Instanton size- dependence of the real spectrum of the Wilson-Dirac operator in 
an instanton background. 

Approximate instantons on the lattice can be obtained either through a cooling 
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procedure [21], or by taking an appropriate lattice transcript of a continuum instan- 
ton field [22, 17]. We will focus on the latter case. In this case numerical studies 
have shown that the real eigenvalues of arc well localised around 0, 2, . . . , 2(i (or 
equivalently, the crossings of zero by eigenvalues of H{m) occur close to these values) 
when the instanton is large at the scale of the lattice spacing, but become delocalised 
as the instanton size is decreased [17]. The standard explanation of this is that in- 
stantons which are small at the scale of the lattice spacing are not slowly varying 
at this scale in the region in which they are localised, so their lattice transcripts are 
"rough" in this region. On the other hand, large instantons are slowly varying, so 
their lattice transcripts are "smooth". The bounds (2.2) can be used to give a more 
precise version of this intuitive explanation as follows. A continuum instanton field 
centred at x^"^ has the form [23] 

where r]'^^ is the 't Hooft symbol, are generators of the SU(2) subgroup and the 
parameter p specifies the size of the instanton. Its curvature is 

F..W--4<- (|^_^,g|,^^). «° (5.2) 

When putting the instanton on the lattice with periodic boundary conditions it is 
important to transform (5.1) to a singular gauge before taking the lattice transcript 
(and the lattice volume must also be sufficiently large that the singular gauge instan- 
ton is close to vanishing at the boundary) [22, 17]. ||F^,^(a;)|| is not affected by this 
though since it is gauge invariant. Prom (5.2) we see that ||F^,^(x)|| diverges at in 
the limit of small instanton size p. Hence for small p the lattice transcripted field gen- 
erally violates the smoothness condition (1.1) since ||1 — ?7(p)|| = \\a'^F^y{x) +0{a^)\\ 
becomes large for plaquettes p close to x^^\ (This is assuming there is no special 
cancellation between a^F^,^(a;) and the O(a^) term; generically there is no reason to 
expect such a cancellation to occur, and in particular when the lattice spacing is small 
a'^Fni,{x) will dominate the 0{a^) term.) Then the locahsation result of §2 on the 

real spectrum of breaks down.^ 

^More precisely, the assumptions under which the localisation was derived break down. This 

19 



On the other hand, from (5.2) we get a bound 




(5.3) 



showing that ||F^,^(a;)|| vanishes uniformly in the hmit of large p. Consequently, for 
large p the smoothness condition (1.1) will generically be satisfied on sufficiently fine 
lattices, thereby guaranteeing localisation of the real spectrum of Du, according to 
the result of §2. 

These considerations can be extended to more general gauge fields describing a 
collection of topologically charged "lumps" (e.g. instanton-anti instanton configura- 
tions, multi-instantons, instanton gases). The topological charge of a lump is given 



If the lump size is small then ||F^j,(x)|| must be large in the lump region in order 

that the magnitude of the integral in (5.4) can be ~ 1. The smaller the lump is, 

the larger ||F^,^(x)|| must be in the lump region. This generically leads to violation 

of the smoothness condition (1.1) as before. On the other hand, if the lump size is 

large ||F^,^(x)|| is not forced to be large in any particular region. Generically we can 

expect ||F^,^(a;)|| to decrease with increasing lump size, and to vanish in the large 

lump limit. Then, by the same argument as before, localisation of the real spectrum 

of Dyj will generically hold in gauge backgrounds describing topological lumps when 

all the lumps are sufficiently large and the lattice is sufficiently fine, 
does not necessarily imply that the localisation result itself must break down, although it is not 
surprising that it should do so. We can turn things around and interpret the numerical results on 
the delocalisation of the real spectrum in small instanton backgrounds as indicating that, in general, 
a smoothness requirement of the form (1.1) is not only sufficient but also a necessary requirement 
for the real spectrum of Dyj to be localised. 



by 




lump 



(5.4) 
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6 Summary 



We have derived general lower bounds on the magnitude of the spectrum of the 
Hermitian Wilson-Dirac operator: 

\H{m)\ > s/l-Cke-\k-m\ for A; = 1, 3, . . . 2d - 1. 

where e is the constraining parameter in the smoothness condition (1.1) (and the 
Wilson parameter is r — 1; the generalisation to arbitrary r > is given in the 
appendix.) Thus we have supplemented the previous bounds for the "physical" case 
k = l [3, 4] with bounds for the "doubler" cases A; = 3, 5, . . . , 2(i — 1. The bounds were 
shown to hold with 

Cfc = c' + 4 

c' = (l + V2)d(d-l)/2 , 4^2<^-^{d-l){d + 2){{k-df -l + d) 

The bounds are rather weak due to the large size of , which is due to the large 
number of terms in the expression (3.19) for x(A;). In practice it is often possible to 
get sharper bounds (i.e. smaller c^) by considering simplified expressions for x(A;). 
For example, in dimension rf=4 we saw how such simplifications lead to bounds with 
c'/ = = 12 and c'l = c'l = 42. In the k = 1 case this is the same as the value 
obtained in [3]. It is plausible that bounds with even smaller c'l can be derived by 
an extension of the arguments of [4] (which gave c'/ = 6) but we did not pursue this. 
For the applications considered in this paper it suffices simply to show that bounds 
of the above form exist, without necessarily finding sharp ones. 

As discussed in §2, the lower bounds on \H{m) \ imply a localisation result on the 
real eigenvalues of the usual Wilson-Dirac operator: the eigenvalues of D^, (in units 
of 1/a) are locahsed around the values 0, 2, 4, . . . , 2(i when e is sufficiently small. (A 
precise formulation of this statement was given in §2.) 

The bounds allow previous results on the overlap Dirac operator to be extended 
from the < m < 2 case to general values of m (m 7^ 0, 2, ... , 2d). This includes 
evaluation of the classical continuum limit of the axial anomaly and index [10, 11], and 
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the results of [3] on locality of the overlap Dirac operator and its smooth dependence 
of the gauge field. The bounds were also seen to imply the existence of topological 
phases for the overlap Dirac operator when one restricts to the space of lattice gauge 
fields satisfying (1.1) with e < 1/cfe for all A;. A complete description of the topological 
phase structure was obtained by combining the bounds with the heuristic result of §4 
on the spectral flow properties of H{m). 

Finally, we pointed out how the bounds can be used to get a more precise un- 
derstanding of why the real spectrum of the Wilson-Dirac operator in an instanton 
background is generally localised around 0, 2, . . . , 2d when the instanton size is large 
but becomes delocahsed when the instanton is small at the scale of the lattice spac- 
ing. (The argument also applies to more general gauge fields describing a collection 
of "topological lumps".) Our argument for delocalisation of the real spectrum in 
small instanton backgrounds involved an assumption, namely that, generically, the 
smoothness condition (1.1) is not only sufficient but also a necessary condition for 
localisation of the spectrum. Numerical studies (e.g. [17]) seem to indicate that this 
is the case, but it would be interesting if it could be proved analytically. This is rele- 
vant for the issue of chiral symmetry breaking in lattice gauge theory since it means 
that the contribution to the density of near-zero eigenvalues of the Dirac operator 
from gauge fields describing small topological lumps is reduced on the lattice. Is this 
reduction an unwanted lattice artifact, or is it a genuinely physical feature revealed by 
lattice regularisation (in the same way that lattice- and other regularisations reveal 
the presence of anomalies which one would not have expected from formal continuum 
considerations)? 
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Appendix 

A Bounds in the case of general Wilson parameter r > 

Using (3.6)-(3.8) a simple calculation gives [5] 

H{m, rf = H{Q, - 2mr^ + r^rn^ (A.l) 

It follows that 

\H{m,r)\ > rm when m<0 (A. 2) 

It is well-known that a lower bound on \H{m) \ is also a lower bound on \H{2d — m)\, 
hence (A. 2) implies 

\H{2d + m,r)\ > rm when m>0 (A.3) 

To see this explicitly, write H{m,r) out according to the definitions (3.1) and (3.8): 

H{m, r,U)=^,[rd-rm + J2 U^'iT+^. - T_^) - r{T+, + T_^))) (A.4) 

It follows that 

H(2d — m,r,U) — —^^(^rd — rm 

+ E Um-T^,) - i-T-,)) - r{{-T^,) + (-T_,))) ) 

(A.5) 

Since (T+^)xj/ = Uf,{x)6a:,y-f, and (T_^)^y = f/^(x-/i)"^4,j/+A the replacement T±^ 
—T±fj_ is equivalent io U ^ —U. Hence (A.5) can be written as 

H{2d-m,r,U) = -H{m,r,-U) (A.6) 
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The operator Rn{U) = 1 — 1(7+/^ + TL^) remains positive under U — > —U, so the 
argument leading to (A. 2) remains vahd under this replacement and we get \H{2d + 
m, r, ?7)| = \H{—m,r,—U)\ > rm for m > as claimed in (A. 3). 

It remains to derive the generalisation of the bounds (2.2) in the general r case. 
In this case the relation (3.10)-(3.13) becomes 

H{m, ry = Y,Sl + r\-m + Y: R^f + E'{r) (A.7) 

where E\r) = Ef^M^'^Y^lSf, , S^] +in^'[S^, C^]) has a bound ||^'(r)|| < c'(r) e. 
A simple generalisation of the argument in [4] shows that this bound is satisfied with 
c'(r) = (1 + rV2)d{d - l)/2. Following [16, 4] we also note that for an eigenvalue 
A(m, r) = {il>{m,r) , H{m,r)il>{m,r)) we have ■£^X{ni,r) — —r{il!{m,r) , 75'0(m, r)) 
and consequently \-£^X{m,r)\ < r, which implies \H{m,r)\ > \H{m',r) \ — r|m' — m|. 

We consider the cases r < 1 and r > 1 separately. In the former case (A.7) 
together with (3.10) gives 

> r^{l + x{m))-c'{r)e (A.8) 

This together with (3.15) gives H{k, r)^ > r^(l — (c^ + d{r)/r'^) e), and consequently, 
setting Cjfc(r) = + c!{r)/r'^ , 

\H[m,r) > r^Jl — Ck{r) e — r\k — m\ , k — 1,3, . . . ,2d — 1 (A. 9) 

In the r > 1 case we rewrite (A.7) as 

> l + x{m) - c'{r)e (A.IO) 
and it follows from (3.15) that 

\H{m,r)\ > ^/l^^d^(r)~e-r\k -m\ , A; = 1, 3, . . . , 2d - 1 (A.ll) 
with Cjk(r) = c'^ + c'{r). 

Note that the bounds (A. 9) for the r < 1 case and (A.ll) for the r > 1 case are 
both weaker than the bound (2.2) for the r = l case. 
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